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ABSTRACT 

We introduce and test an expression for calculating the variance of a physical field in 
three dimensions using only information contained in the two-dimensional projection 
of the field. The method is general but assumes statistical isotropy. To test the method 
we apply it to numerical simulations of hydrodynamic and magnetohydrodynamic tur- 
bulence in molecular clouds, and demonstrate that it can recover the 3D normalised 
density variance with 10% accuracy if the assumption of isotropy is valid. We show 
that the assumption of isotropy breaks down at low sonic Mach number if the tur- 
bulence is sub-Alfvenic. Theoretical predictions suggest that the 3D density variance 
should increase proportionally to the square of the Mach number of the turbulence. 
Application of our method will allow this prediction to be tested observationally and 
therefore constrain a large body of analytic models of star formation that rely on it. 

Key words: ISMxlouds - ISM: kinematics and dynamics - turbulence - magneto- 
hydrodynamics - methods: statistical. 



1 INTRODUCTION 

It is a fundamental problem in astrophysics that we typ- 
ically only have access to two-dimensional (2D) physical 
fields that have been integrated or averaged over the line- 
of-sight, while the physical fields of interest are intrinsically 
three-dimensional (3D). Constraining 3D properties is very 
difficult, and is limited to structures with a fairly simple ge- 
ometry (e.g. Lucy 1974 - see Reblinsky 2000 for an applica- 
tion to galaxy clusters). Molecular clouds have very complex 
structure that is not suitable for direct inversion from 2D to 
3D, but there exist methods for inferring statistical infor- 
mation on the scaling behaviour of the 3D structure of the 
density and velocity fields (e.g. Stutzki et al 1998). Molec- 
ular cloud evolution is driven in large part by the action 
of turbulence, and recently, much interest has been directed 
towards the action of turbulence in shaping the density field 
in molecular clouds, in particular towards the theoretically 
predicted increase in the variance of the density field with 
the Mach number of the turbulence (Padoan, Nordlund, & 
Jones 1997; Passot & Vazquez-Semadeni 1998). This pre- 
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diction plays a key role in analytic models of star forma- 
tion (Padoan & Nordlund 2002; Krumholz & McKee 2005; 
Elmegreen 2008; Hennebelle & Chabrier 2008; Padoan & 
Nordlund 2009; Hennebelle & Chabrier 2009). 

If the density, p, is expressed in units of the mean den- 
sity, po, the theoretical predictions for the relationship be- 
tween the density variance and root mean square (rms) Mach 
number, M, may be simply written as: 

= 6^M^ (f) 

where & is a constant of proportionality. In the case of 
isothermal gas, which is applicable to a good approxima- 
tion to molecular clouds, the probability density function 
(PDF) of the density field is thought to be lognormal in 
form (Vazquez-Semadeni 1994). Using a lognormal form for 
the density PDF, Padoan & Nordlund (2002) derive a re- 
lation between the variance in the logarithm of the density 
field, o"hi(p/p(,), and the rms Mach number, M: 

'^ln(p/P0) + (2) 

which is equivalent to equation ([T}. A range of predictions 
for b have been proposed, which were recently synthesized 
into a unified model by Federrath et al (2008) who propose 
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that b — 1/3 for solenoidal (divergence- free) forcing and 
b — 1 for compressive (curl- free) forcing in 3D. 

Given the importance of equation for analytic mod- 
els of star formation, it is essential to test it with observa- 
tional data. Initial observational tests by Goodman, Pineda, 
& Schnee (2009) did not find any obvious support for the 
theoretical predictions. Federrath et al (2009) have sug- 
gested a reason for this lack of agreement, citing variations 
in b caused by different turbulent forcing mechanisms. Per- 
haps more importantly, these observations relied on measur- 
ing the variance in the projected 2D column density rather 
than in the 3D density field that appears in the theoretical 
predictions. From observations, we do not have access to the 
3D density and velocity fields to directly test this prediction, 
and evaluate its applicability to molecular clouds and to the 
predictive star formation models. 

In this paper, we derive and test an expression that re- 
lates the observable variance in the 2D column density field 
to the true 3D variance of the density field. The method is 
completely general, albeit with the assumption of statisti- 
cal isotropy, and may also be applied to the projection of 
other physical fields. In Section 2, we present the analytical 
expressions needed to convert the 2D variance into the 3D 
variance. In Section 3, we demonstrate the application of 
the method to numerical simulations of hydrodynamic and 
magnetohydrodynamic turbulence. Our summary is given in 
Section 4. 



2 THE STATISTICS OF PROJECTION FROM 
3D TO 2D 

2.1 Development of the Method 

Given a physical field in 3D, we wish to examine how the 
statistical properties of a 2D projected field are related to the 
intrinsic properties. By doing this, we can then infer intrinsic 
quantities from the measureable quantities. We focus here 
on the variance, measured in 3D and in projected 2D. In our 
initial derivation, we consider a 3D field, Fi(x,y, z), which 
is averaged along the z axis to produce a projected field 
F2{x,y), via: 



1 r^/^ 

F2{x,y) = - dz Fi{x,y,z) 

^ J-L/2 



(3) 



where L is the physical size of the region, assumed cubical, 
that contains Fz. 

We will make use of the 3D Fourier series of ^3 , obtained 
over the interval [—L/2,-\-L/2], which is: 



L/2 fL/2 fL/2 



L/2 J-L/2 J-L/2 



d^r F'iir) e 



(4) 



where r = (xx,yy,zz) and the spatial frequencies are fe = 
{kxX,kyy,kzz), for integer kx, ky, and kz- The field F3 can 
then be written as: 



..00 00 00 

^3^ = 13 E E E Mk.,ky,kz 



27Tik-r / L 



kj- — — 00 fci, — — 00 - — — 00 



(5) 



Inserting equation ((5]) into equation ([3]), 

F2{x,y) = 
-, rL/2 00 00 00 

JlJ dz E E A(fc.,fc.,fc.)e'^"=-'-/^(6) 

and computing the z integral first, we find that: 



00 00 



F2{x,y)^^ J2 E Mk.,ky,kz^O)e 



1-Kik-r j L 



krr — — 00 fc,, — — 00 



Since 

L J-L/2 
1 /■^/2 



dz e 



27vik^z / L 



= 1 for kz = 0, 



dz e 



2'KikzZ / L 



L/2 



for kz ^ 0. 



(7) 



(8) 



(9) 



Since the 2D Fourier series, over the interval [—L/2, -fL/2], 
of F2 is: 



(•L/2 pL/2 

F2{k2) = I / d'ra F2(r2) e 

-L/2 J-L/2 



2iTik2-T2/L (10) 

where r2 — {xx,yy) and = {k^x, kyij), and: 

(11) 



F2{r2) = ^ ^ ^ F2(k^,ky)e'- 

fca; — — 00 ky — — 00 



k2-r2/L 



we see that, comparing equation (0 and equation (lll|l : 

F2{kx,ky)^ ^F3{kx,ky,kz^0). (12) 

In other words: the 2D Fourier series of F2 is propor- 
tional to the kz — cut through L3. Previous studies have 
made use of this result (e.g. Stutzki et al. 1998; Lazar- 
ian & Pogosyan 2000; Brunt & Mac Low 2004; Miville- 
Deschenes & Martin (2007)). As the power spectrum, P3(fc), 
is the squared modulus of the Fourier transform {Pz{k) — 
F3{k)F3 (fe)), 2D power spectra of projected fields have 
been used to infer the 3D power spectrum, under the as- 
sumption of isotropy, i.e. that the kz = cut through the 
power spectrum is statistically representive of the full power 
spectrum. In our analysis below, we will also make use of the 
power spectrum and the assumption of isotropy. 

Now we compute the variances in the 3D and 2D fields. 
From the mean value of F^ : 



L3 



L/2 fL/2 rL/2 



d^rFz(r) 



-L/2 J -L/2 J -L/2 

and the mean square value of L3: 

fL/2 fL/2 fL/2 



d'rF,'{r), 



-L/2 J-L/2 J-L/2 

we find the variance of F3 as: 



= {F^} - {Fsy 



(13) 



(14) 



(15) 



We also make use of the Fourier transform of L3 , by noting 
that: 



(Fa) = —i?3 (0,0,0) 



(16) 
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and, through Parseval's Theorem: 

^ oo oo oc 

= ^ E E E ^3F3*, 

fca^ — — OO fcy — — cxD k~— — oo 

allowing us to rewrite equation (|15p as: 



(17) 



2 
0"3 



^ E E E ^3i?3* U7?/(o,o,o) 

\ \ fcx — — oo fey — — oofc^ — — oo / / 



(18) 



By similar analysis, we find the variance of F2 is: 



2 

0-2 



IT ( ( E E ^2i?2* ) -i?2'(0,0) j , (19) 



or, applying equation ([12 



2 1 
-2 = ;^x 



00 00 



E E ^3(fc. =0)F3*(fc. = 0) -#3^(0,0,0) . 

y yfcj;— — 00 fcy— — 00 y y 

(20) 

The ratio, R, of the variance of F2 to the variance of F3 is 
therefore; 

77 = R = 

(00 00 \ 

E E i?3(fc. =0)i?3*(fc. = 0) -i?3'(0,0,0) 
fca;— — 00 fc^y— — 00 / 



E E E ^3^3* - 



F3'(0,0,0) 



1 — — 00 ku = — oo k, — - 



(21) 

In practice, the measureable physical fields, either 
through observations or numerical simulations, will consist 
of a discrete number of measurements at a fixed set of grid 
points. We define the scale ratio, A, as the ratio of the image 
(or cube) size to the pixel size. Taking the grid to be of size 
AxAxAorAxA pixels in 3D and 2D respectively, the Fourier 
transforms are carried out at a discrete set of spatial frequen- 
cies, k = -A/2 + 1, -A/2 + 2, -2, -1, 0, 1, 2..., A/2-1, A/2 
along each axis. The quantity -F3-F3 is the spectral power, 
Pa, which can be observationally derived (up to an unim- 
portant constant of proportionality) as follows. If we had a 
projected mean 2D field, F2, and wanted to use this to infer 
the 3D variance, then we would compute its power spectrum, 
P2{kx,ky), and from this produce an azimuthally-averaged 
power spectrum P2(k), which depends only on the modulus 
of the spatial frequency, k. The key idea behind the method 
is that we can, assuming isotropy, take: 

Paik) (X P2ik), (22) 

to obtain: 

/ A/2 A/2 \ 

E E P2ik)]-P2{0) 

\k:c = -\/2+l ky = -\/2+l I 



R 



' A/2 A/2 A/2 \ 

E E E pm]~pm 

, fc^ = -A/2+l fc„ = -A/2+l fc^ = -A/2+l / 



(23) 

The 3D variance can then be calculated as (jf — cr| / R. Note 
however that since A is necessarily finite, the observed 2D 
variance and the calculated 3D variance are lower limits to 
the true variances that would be obtained in the limit A — > 
00. This is discussed further below. 

A more compact form of equation (|23|l can be obtained 
by defining: 



E 

fc5^0 

and; 



(A/2 A/2 ' 

E E 
fex = -A/2+l fcy = -A/2+l , 



-P2(0) (24) 



E ^2(fc) 

fc#0 
/ A/2 

E 



A/2 A/2 > 

E E ^^(^) 



P2(0) 



-A/2+1 fcy = -A/2+l fc^ = -A/2+l 



SO that: 



R 



E 

E ^2(fe) 



(25) 



(26) 



For this method to work, it is essential that the pro- 
jected field, F2, is the line-of-sight average of the 3D field, -F3 
(see equation (O). In many instances, the projected 2D field 
is the line-of-sight intregral of the 3D field (e.g. column den- 
sity versus density). A simple solution for integrated fields 
that ensures the above requirements are satisfied is to ex- 
press the field F2 in normalised units ~ i.e. by dividing F2 by 
its mean value. In this case, the variance of the normalised 
3D density field can be calculated from the variance of the 
normalised 2D column density field, as discussed below in 
Section 2.6. 



2.2 Approximations 

The power spectra of many fields of interest are power-law 
in form; P3 cx k~" . If the power spectrum is steep (large 
q) then the variance is sensitively dependent on the power 
at low spatial frequencies. Because of the quantisation in 
equation (I26p and the small amount of information available 
at low k, the above procedure may be inaccurate for large 
Q. Brunt & Mac Low (2004) examined 3D and projected 2D 
standard deviations in the normalised density field (p/po) 
and normalised column density field (N/Nq) obtained from 
numerical simulations, with A — 128. The (column) density 
power spectra were reasonably well fitted by Pi{k) oc k'"^, 
albeit with some curvature. They found that un/Nq/'^p/po ~ 
0.34 ± 0.04. Evaluating equation ((26j, using A = 128 and 
Psik) oc k~^, to compute the ratio of standard deviations 
(i.e. ^/R), we predict 0-2/0-3 = 0.39, which is in acceptable 
agreement with the Brunt & Mac Low (2004) experimental 
result, given the spectral curvature. 

With a power law form for the power spectrum, it is 
tempting to approximate equation (|26p with simple integral- 
based expressions for R. If the spectral slope is a we could 
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write: 

2tt r^/^ dk k k-"' 

R= ^7^ , 27 

4nJ^' dk fc2 fc-" 

which is easily solved for a specified a. However, this is not a 
terribly good approximation as the integrals have been com- 
puted over circular and spherical regions of fc-space rather 
than the square and cubical regions which contain the fields. 
At low a this is not sufficiently accurate, and in general we 
recommend that the direct summation method presented 
in equation (|26p be used. As a trivial example, for a cube 
of Gaussian noise (q = 0), equation H26|) correctly predicts 
R = (X^ ~ 1)/(A^ - 1), which tends to 1/A for large A. In 
the same limit, equation (|27|) incorrectly predicts R — 3/2A. 
Equation (|27|) is useful however in gaining some intuitive 
understanding of how R depends on the form of the power 
spectrum: R decreases if proportionally more power resides 
at large k because the denominator is weighted by k^ and 
the numerator by k. The physical effect underlying this is 
that, if the variance is mostly at high k (small spatial scales), 
line-of-sight averaging suppresses more of the power than if 
the variance is mostly at low k (large spatial scales). Equa- 
tion (|27|) is also useful for exploring scale-dependence effects 
on the calculated variances. Note that for 2 < a 3, the 2D 
variance converges as A — y oo but the 3D variance diverges! 

2.3 Application to Observational Data 

For fields acquired through observations, it is necessary to 
account for the telescope's beam response (point spread 
function) in the calculation of variances. The observational 
version of equation (I26p is: 



E P2ik)BHk) 

where A is the number of pixels along each axis, -B^(fc) — 
B{k)B* (k) is the square of the Fourier space representa- 
tion of the telescope beam pattern, and P2{k) is the power 
spectrum of F2 in the absence of beam-smoothing and 
instrumental noise. (Note that the observable quantity is 
P2{k)B^{k) + P2N{k) where P2jv(fc) is the noise power spec- 
trum.) Equation (I28|l includes accounting for the fact that 
the observed variance in the 2D field has been suppressed 
by the smoothing imposed by the telescope beam, and this 
must be taken into account. In equation (|28p . we have also 
applied the beam pattern to the denominator (represent- 
ing the 3D variance). This procedure may at first appear 
an odd choice, but it is motivated by the requirement that 
we limit our knowledge of the 3D variance to the resolution 
provided by the data. In effect, equation (|28|l is an attempt 
to construct a 3D beam that samples the 3D density field 
and affords the same effective linear resolution as that pro- 
vided by the telescope beam in projected 2D. The Fourier 
series F3 extends to fc = cxa, or at least to a (potentially 
very large) spatial frequency, kc, beyond which there are no 
variations in the field. In other words, there will be fiuctua- 
tions in the projected 2D field (and therefore the 3D field) 
that we are not directly sensitive to, as they lie below our 



resolution limit. These, naturally, are additional sources of 
variance that we cannot measure. 

The full variance of the field Fs should in principle be 
derived by removing (from the denominator) the effect of the 
beam pattern, and summing over the spatial frequency range 
1 ^ k !^ kc. In general, we will have little or no information 
about kc, unless of course a distinct break in the power spec- 
trum is observed. A possible, theoretically-motivated choice 
is that kc is the spatial frequency corresponding to the sonic 
scale (Vazquez-Semadeni et al 2003; Federrath et al 2009). 
We may try to account for this, by assuming a form for P3 
in the inaccessible range kc > k > A/2, by extending an 
observed power-law for example. This is quite dangerous, 
particularly for a ^ 3 for which the summation diverges at 
large spatial frequencies. The case a — 3 provides a con- 
venient reference point, at which R is underestimated by a 
factor ~ In (A/2)/ In (fee) - a logarithmic divergence. With- 
out a direct constraint on kc it is more reasonable to restrict 
our knowledge of the 3D variance to the spatial dynamic 
range afforded by the observations, and utilize equation ([28} 
as stated, while recognizing that the 3D variance must be a 
lower limit to the true 3D variance. For numerical simula- 
tions, this is not an issue, as we have access to the full range 
of spatial frequencies at which structure is present. 

2.4 Accounting for Non-Periodic Fields 

The Fourier series utilised above implicitly assume that the 
field is periodic. While this is true of some numerical simu- 
lations, it is obviously not true of real physical fields. If the 
field F2 has a significant amplitude at the field boundaries, 
then elevated power is produced in the 2D power spectrum 
at some spatial frequencies, typically along the fc^, = and 
ky — axes, and this is not consistent with an isotropic form 
for the 3D power spectrum. Firstly, we note that a signif- 
cant amplitude of F2 at the field edges is almost certainly 
not consistent with the assumption of isotropy, as it implies 
that F2 has been extracted from a larger field with signifi- 
cant power on larger scales. In this case, restriction of the 
field in the projected x and y directions is imposed but there 
has been no corresponding restriction on the extent of the 
field in the z direction. Sensible definition of the field F2 is 
therefore recommended, with, ideally, F2 falling to near zero 
amplitude at the field boundaries. 

It is possible to ameliorate the effect of edge disconti- 
nuities by edge-tapering the field (e.g. Brunt & Mac Low 
2004) . In this case we recommend that A is sufficiently large 
that the tapering causes a negligible modification of the 2D 
variance. A good solution is, in addition to tapering, to pad 
the field with zero values. If an observed field with scale ratio 
A is zero-padded out to a larger size Ap, then it is straight- 
foward to show that the 3D variance should be calculated 
via: 

= + {F2fW - {F2f)/Rr, + {F2?) - {F2)\ (29) 

where r] = Ap/A and Rp is the 2D-to-3D variance ratio cal- 
culated from the power spectrum of the padded field. The 
quantities {F2), o"!, and 0-3 apply to the unpadded field. 

Above we have assumed a square (cubical) box that con- 
tains the field F2 {F3). In practice, we may wish to apply 
this method to fields that are not exactly square, although 
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this then obviously raises questions about the assumption 
of isotropy. For a field of size A^; x Ai,, we recommend zero- 
padding to produce a square field of size Xp^ x Xpy , and using 
equation with rj — [XpxXpy / \x\y)^ ■ Note that this as- 
sumes that the line-of-sight extent of F3 is Az = (\x\y)^ ■ 
Clearly, one should try to ensure that A^, « \y to respect 
the assumption of isotropy. 

2.5 Summary of the Method 

To summarize the method: to estimate the 3D variance, cri, 
the following procedure should be followed: 

(1) measure the variance, al, of the normalised projected 
field F2, 

(2) measure the power spectrum, P2{k), of F2, and 
assume, through isotropy, that P-i(k) oc P2(k), 

(3) using P2{k), compute the 2D-to-3D variance ratio, 
R, via equation (|26|) or equation (|28p . 

(4) compute (t|, accounting for any zero-padding used to 
compute the power spectrum (equation (|29|) ). 



2.6 Density Fields 

So far, other than isotropy, we have not imposed any partic- 
ular structure to the field ^3 (e.g. by choosing a form for the 
power spectrum and PDF). It is well established that pro- 
jected column density power spectra of molecular clouds are 
power-law in form (e.g. Stutzki et al 1998; Bensch, Stutzki, 
& Ossenkopf 2001), implying also that the 3D power spectra 
are also of power-law form. Typically, the 2D power spec- 
trum is of the form P2 oc where a ~ 3. Density fields 
are always positive, and it is convenient to express both the 
density, p, and column density, N, in normalised units, by 
dividing by their respective mean values, po, and A'o, re- 
spectively. In this way, they will conform to the properties 
of the fields ^3 and F2 discussed above (i.e. that F2 is the 
line of sight average of F3). Otherwise, the column density 
will be scaled by the physical length, L, as it is the integral 
of p along the line-of-sight, rather than the average. In most 
cases, our 2D observations will be of column density - e.g. 
through an optically thin spectral line, or through extinc- 
tion mapping. The physical line-of-sight length, L, needed 
to convert column density to projected mean density, may 
be unknown. While L may be inferred through the assump- 
tion of isotropy, if the distance is known, it is better to use 
the normalised density, p/po, and column density, N/Nq. 

Using the observed variance, c^/jy^, in the normalised 
column density field, and the angular-averaged power spec- 
trum, -Pjv/jVo(^)i we can then use equation (|28p to obtain 
R. With appropriate treatment of any zero-padding (equa- 
tion dm) with (Fa) = 1, (t| = a%/j^^, and erf = cr^/pj, the 
variance in normalised density, o-p/p^, can be calculated. 

If a lognormal PDF in 3D is assumed, we can then derive 

'^ln(p/po) Vi^'- 

'^fn{p/po) = In (1 + C^p/po)- (30) 



Note that it is not necessary to assume a lognormal form 
for the column density PDF, although this indeed may be 
true (Ostriker, Stone, & Gammie 2001; Vazquez-Semadeni 
& Garcia 2001). 

Considering that Cp/pg may be subject (for spectral 
slope a ~ 3) to a logarithmic divergence because of unre- 
solved density structure, as discussed in the previous section, 
equation (|30|l has the fortunate property that the computa- 
tion of o"hi(p/p(,) suppresses this divergence to a log(log) di- 
vergence, so that o-fn(p/po) can be quite well estimated even 
from observations at finite resolution, provided the spatial 
dynamic range is sufficiently large to measure c^r/jvo 
curately, and the field is of good linear spatial resolution 
such that A/2 approaches reasonable expectations for k^. 
We recommend that estimates of the unresolved variance 
be made, using equation (|27p as a guide for suitable values 
of kc (e.g. the spatial frequency corresponding to the sonic 
scale, or even kc — > 00) to better assess the utility of the 
observational measurement (Brunt 2009). We note also that 
the above method for estimating o'p^p^ does not require the 
cloud distance to be known. However, observations at fixed 
angular resolution applied to clouds at different distances 
will probe different physical scales, so the method is not 
entirely distance-independent. 

The above prescription can, in principle, be applied to 
other positive-definite fields, such as temperature. Estima- 
tion of mean projected temperature fields is usually done by 
quite a different procedure than that used for column den- 
sity fields, however. Typically, we are able to derive a line-of- 
sight average temperature for each pixel, either through fiux 
ratios at far infrared wavelengths, or by excitation analysis 
of millimetre-wave spectral lines, for example. These mea- 
surements are often not straightforward spatial averages, as 
assumed by our method, but may be (e.g.) density-weighted 
averages instead. Some caution must be applied to the treat- 
ment of such fields. 



2.7 Velocity Fields 

Application of the method to projected mean velocity fields 
is also possible. Obviously, since a velocity field is not 
positive-definite, we do not normalize the field as was done 
for the density field. In principle, we do have access to the 
mean velocity field, which can be obtained through an op- 
tically thin spectral line. However, such a field is density- 
weighted, rather than a direct spatial average as assumed 
by equation Q, and significant problems in estimating the 
power spectrum (Brunt & Mac Low 2004), and presumably 
therefore the 2D variance, can arise for supersonic turbu- 
lence. There will be lines of sight with insufficient signal- 
to-noise to calculate a mean velocity. Therefore the variance 
must be calculated for the detected lines of sight, with the as- 
sumption that these reliably represent the entire field. Note 
that the problem in estimating the mean of the field (as 
required in the case of column density / density analysis) 
does not apply here. We also have access to only one com- 
ponent of the velocity field (the line-of-sight component v^) 
and we must assume not only isotropy in this, but also ap- 
ply the isotropic assumption to the entirely inaccessible Vx 
and Vy components. With these provisos, the method would 
proceed as outlined at the end of Section 2.1. 



6 C. M. Brunt, C. Federrath, & Price, D. J. 




1 10 




Figure 1. Plot of the estimated 3D density variance, a-^/i?, ver- 
sus tlie actual 3D density variance, o-g, obtained from numerical 
simulations of hydrodynamic turbulence. Each plotted point rep- 
resents the mean value of /R obtained from analysis of column 
density fields projected along the three cardinal directions, while 
the error bars represent the standard deviation in the values of 
(Tj/R for the three directions. 
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Figure 2. Plot of the estimated 3D density variance, (T|/iJ, ver- 
sus the actual 3D density variance, o-g, obtained from numerical 
simulations of magnetohydrodynamic turbulence. Each plotted 
point represents the mean value of cr|/R obtained from analysis 
of column density fields projected along the three cardinal direc- 
tions, while the error bars represent the standard deviation in the 
values of U2/R for the three directions. The solid line is the line 
of equality. 



There is a rather more straightforward method of es- 
timating the line-of-sight velocity dispersion, (t„. , which is 
simply to measure the dispersion of the mean line profile 
obtained through imaging observations of an optically-thin 
spectral line. (This is also a density-weighted measure, how- 
ever.) If we have an estimate of the cloud temperature, and 
an estimate of the mean molecular weight, we can then de- 
rive the Mach number to construct the right-hand-side of 
equation ([T|. 



3 APPLICATION TO NUMERICAL 
SIMULATIONS 

As an initial test of the method, we now apply it to numeri- 
cal simulations of hydrodynamic and magnetohydrodynamic 
turbulence. Given that the above theory is derived assuming 
isotropic fields, we investigate the effect of anisotropy using 
magnetohydrodynamic simulations. In all the studies below, 
we use normalised column densities (N/No) and normalised 
densities (p/po) so that (t| = o"Ar/]VQ ^nd 
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3.1 Hydrodynamic Turbulence 

First, we analyse hydrodynamic simulations at a range of 
rms Mach numbers: 1.25, 2.5, 3.5, 5, 7, 10, and 20. Multiple 
snapshots are used at each Mach number setting, separated 
by at least a crossing time. The hydrodynamic simulations 
were run with the phantom Smoothed Particle Hydrody- 
namics code, with turbulence driven artificially over 5 cross- 
ing times using large-scale solenoidal Fourier driving (Price 
& Federrath 2009; Federrath et al 2009; Brunt, Heyer, & 
Mac Low 2009; Brunt 2003). For the analysis in this paper 
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Figure 3. Plot of the standard deviation in 1J2/R expressed as 
a fraction of the moan (t|/_R, versus the Alfvenic Mach number. 
Ma, for the MHD simulations. 



we performed a suite of low resolution 128^ particle calcula- 
tions, interpolated to 256^ grids. Since the key idea here is 
the reconstruction of the density variance, we find that even 
low resolution calculations are sufficient as demonstrators 
of the method. For each simulation, we take the 3D density 
field, of mean unity, and create three normalised projected 
fields by averaging along each spatial axis. We compute ct| 
in 3D and a value of af for each projected 2D field. 

The power spectrum for each column density field is 
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Figure 4. Relations between o"|/iJ and the true for each of the three axes. Runs with ^ 1 are shown as dots; runs with AI^ > 1 
are shown as open squares. The solid line is the line of equality and the dotted lines mark underestimation and overestimation of Cg by 
a factor of 2. 



then measured (specifically, we compute the average power 
in bins of unit width in k.) Following this, equation (|26p is 
used to compute a value of R for each projected field. To do 
this, we interpolate the binned column density power spec- 
trum to the appropriate power for each k resulting from the 
nested sums. An extrapolation is required (using a power 
law fit to the power spectrum in the high k region) for 
\/2A/2 < fc ^ \/^X/2 in the denominator. In the above pro- 
cedure, we emphasize that only information present in the 
column density field is used in the calculation of R. Values 
of R computed for these fields vary between about 0.03 and 
0.15. Note that in general, R will depend on the form of the 
power spectrum and the scale ratio of the field. 

In figure [T] we plot the estimated 3D density variance, 
CT2/R, versus the actual 3D density variance, erf. For each 
snapshot, we represent the measurements as a mean (plotted 
point) and standard deviation (error bar) obtained from the 
three different projections. The method can predict the 3D 
density variance to about 10% accuracy. 

3.2 Magnetohydrodynamic Turbulence 

We now analyse magnetohydrodynamic (MHD) simulations 
of turbulence. The MHD simulations were run with the grid- 
based code FLASH (Fryxell et al. 2000). We used a new ap- 
proximate Riemann solver for ideal-MHD (Bouchut, Klin- 
genberg, & Waagan 2007, 2009), which preserves positive 
states in highly supersonic MHD turbulence. This solver 
was recently developed for FLASH by Waagan (2009). The 
corresponding scheme for preserving positive states in hy- 
drodynamical studies has been successfully tested and ap- 
phed in Klingenberg, Schmidt & Waagan (2007). We used 
the same solenoidal forcing scheme as used in Federrath et 
al. (2009) and Price & Federrath (2009). All MHD models 
were evolved on a fixed grid with 256"^ grid zones. 

In addition to varying the Mach number, in the MHD 
simulations we now also vary the ratio of thermal to mag- 
netic pressure, /? = P/{y^). For large magnetic field 
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Figure 5. Plot of the anisotropy indicator, o"|^/o"2||i versus the 
Alfvenic Mach number. Ma, for the MHD simulations. For ref- 
erence, the horizontal dashed line marks the expected value of 
o-jj^/ctjii for an isotropic field; the vertical dashed line marks 
Ma = 1, below which anisotropy is notable. 



strengths (lower /3) significant anisotropy can arise in the 
density and velocity fields (Mac Low 1999; Vestuto, Os- 
triker, & Stone 2003; Heyer et al. 2008), and thus the as- 
sumptions of our method could break down. Using the MHD 
density fields, we repeat the same procedure described above 
to calculate R, and hence the estimated 3D variance, erf /R; 
the comparison of estimated to true 3D variance is shown 
in Figure (2] We find that the method is still accurate in 
predicting the mean (plotted points) but the variation be- 
tween the different directions (error bars) is now signifi- 
cantly larger for some parameterizations. Investigation of 
this shows that significant variations in 02/ R between the 
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Figure 6. Representative 2D slices through a sample of density fields created by sub-Alfvenic turbulence. The mean magnetic field 
direction is illustrated by an arrow, labelled Bq. The values of the sonic Mach number, M, the plasma beta, /3, and the Alfvenic Mach 
number. Ma are given above each panel. 




Figure 7. Power spectra of the 2D density field slices from Figure |6] Greyscales are represented in logarithmic units, with contours, 
each separated by a decade in power, overlayed; the contours been smoothed for clarity. 



different projection axes are seen if the Alfvenic Mach num- 
ber, Ma = M\J (/3/2), is smaller than ~ unity. To show this, 
in Figure [3] we plot the fractional error in the estimated 3D 
density variance (i.e the standard deviation \mj\lR divided 
by the mean o^jR obtained from the three different pro- 
jection directions) versus Ma- The fractional error in the 
estimated 3D density variance is around 10% for Ma > 1, 
while it increases significantly for Ma < 1- 

As there is now notable variation in a'i/R between the 
different directions, it is worth plotting the individual o^jR 
versus a\ relations obtained from each axis. These relations 
are shown in Figure |4l where we have distinguished the plot- 
ted points with different symbols for Ma ^ 1 and Ma > 1- 
The 3D variance in sub-Alfvenic runs is overestimated (un- 
derestimated) from 2D fields produced by averaging over 
an axis parallel (perpendicular) to the mean magnetic field 
direction, by as much as a factor of ~ 2. 

The variations in cri /R are caused by anisotropic struc- 



ture in the density field that is produced when Ma < 1. As a 
simple measure of this anisotropy, we define an "anisotropy 
indicator" using the 2D fields that are produced by averag- 
ing over the three cardinal axes. The anisotropy indicator is 
defined as o-2x/o-||| = o\y) I2a\^, where crL, crfy, and 

are the variances in the 2D fields produced by averaging 
over the x, y, and z axes respectively. The measured values 
of o"2x/o"2|| s^re plotted in Figure [S] versus Ma- For Ma > 1, 
'^2±/'^2\\ is ~ unity, indicating roughly equal variances in 
each of the projected 2D fields. For Ma < 1 the variance in 
the 2D field averaged over the z direction (parallel to the 
mean B field) is larger, by as much as a factor of 2, than 
that of the 2D fields produced by averaging over the x ox y 
directions (perpendicular to the mean B field). 

Interestingly, the degree of anisotropy appears to be re- 
lated to the sonic Mach number, M, as is evident in both 
Figure[3]and Figure[5] Indeed, for M > 10, the estimation of 
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Figure 8. Column density plots (logarithmic greyscale) of the non-periodic SPH simulation. 



the 3D variance is not noticeably worse for sub-Alfvenic tur- 
bulence than for super- Alfvenic turbulence, and the largest 
deviations from isotropic behaviour are found in the M = 2 
runs. A direct visualization of the anisotropies present in the 
density fields is given in Figure|6l which shows representative 
2D slices through a sample of density fields produced by sub- 
Alfvenic turbulence. At low sonic Mach number (M = 2) 
the density field is characterized by many long filamentary 
structures oriented parallel to the mean magnetic field direc- 
tion. As the sonic Mach number is increased, this directional 
order in the density field is systematically reduced, until at 
M = 10 it is barely noticeable. The power spectra of these 
2D density field slices are shown in Figure [7] These con- 
firm our visual assessment of the density field slices. The 
high frequency spectral power along the magnetic field axis 
is notably lower in the low sonic Mach number calculations. 

The tests on the MHD simulations provide a valuable 
baseline for establishing the physical regimes in which our 
method can be applied. In most circumstances encountered 
in molecular clouds with sizes of a few parsecs, sonic Mach 
numbers are large enough to enforce sufficient isotropy that 
the method can work with around 10% accuracy, even for 
sub-Alfvenic turbulence. However, these are rather ideal- 
ized conditions, as the models do not include other rele- 
vant physics such as self-gravity. Gravitational collapse in 
the presence of a magnetic field is likely to be oriented along 
the magnetic field lines if the field is sufficiently strong. Note 
that the combination of magnetic fields and gravity, in this 
instance, can induce anisotropic structure in the orthogo- 
nal direction to that seen in the sub-Alfvenic simulations 
described above. 



3.3 Non- Periodic Fields 

In the foregoing analysis, we have used periodic fields for 
which there are no edge discontinuities and consequently no 
zero-padding is required. A more physically realistic scenario 
is provided by SPH simulations that do not employ periodic 
boundary conditions. For our analysis, we take density fields 
produced by the hydrodynamic simulations of Price & Bate 
(2009). Column density fields extracted at three difi'erent 
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Figure 9. Comparison of calculated to true 3D variances for the 
non-periodic SPH density fields. 



times (at t — 0; after one free fall time; and at the end of 
the calculation where t — 1.274 tff) are shown in Figure [S] 
(The column density fields are represented as logarithmic 
greyscales as the variances are extremely large.) We zero- 
pad these fields to Xp = 512 from the original A = 256, then 
compute the power spectra and from this derive Rp. We then 
calculate the 3D variance by use of equation (|29p . 

In Figure]^ we compare the 3D variances derived from 
our method with the actual 3D variances measured directly 
from the density fields. We find that the method is com- 
parably accurate for these non-periodic fields, albeit with 
a rather large variation (~ 40%) between different projec- 
tions in the latest snapshot. In this field, the variance is 
very large indeed (7.7 x 10^) and is dominated by a small 
number of very large density values caused by localised grav- 
itational collapse. The variances in the initial density field 
(a spherical blob), in the t = tff density field, and in the 
t = 1.274 tff density field are recovered by the mean cal- 
culated variances to within 0.3%, 20% and 20% accuracy 
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respectively. To evaluate the success of these measurements, 
we calculate an anisotropy indicator as was done in Sec- 
tion 3.2, by taking the ratio of the variance in one pro- 
jected field to the mean of the variances in the other two 
projected fields. Since the simulations contain no magnetic 
fields, there is no preferred axis. Therefore we combine the 
three 2D variance measurements in a way that minimises the 
anisotropy indicator (lower values represent greater degrees 
of anisotropy). For the three snapshots, in increasing-time- 
order, we find anisotropy indicators of 1.00, 0.78, and 0.54. 
For the two evolved snapshots, the density fields are compa- 
rably anisotropic to the low Ma, low M density fields of the 
previous section, but with the anisotropy now arising from 
the initial velocity perturbations rather than from the effect 
of a magnetic field. The accuracy of the method for non- 
periodic boundary conditions is therefore comparable to the 
periodic case when compared at the same level of intrinsic 
anisotropy. 

The extremely high variances seen in the evolved 
density fields are a result of small-scale collapse in this 
gravitationally-bound cloud, and far exceed variances ex- 
pected by equation IT} with b ~ 0.5. It will be interesting in 
future to see how equation ([!]) is modified in strongly-self- 
gravitating clouds, and whether observations can quantify 
this. 



4 CONCLUSIONS 

We have introduced and tested a simple method for mea- 
suring the 3D density variance in molecular clouds, using 
only information present in the projected column density 
field. The method relies on the assumption of isotropy and 
uses the measured column density power spectrum in con- 
junction with Parseval's theorem to calculate a correction 
factor, R, that scales the observed normalised column den- 
sity variance to the intrinsic three-dimensional normalised 
density variance. The method is sufficiently general to be 
applied to any isotropic field. 

For density fields produced in supersonic hydrodynamic 
and magnetohydrodynamic turbulence, the method is accu- 
rate to about 10% provided that the assumption of isotropy 
is valid. For turbulent density fields, in practice this requires 
that the turbulence motions are super- Alfvenic (Ma < 1), 
though even in the sub-Alfvenic regime we are able to re- 
cover the 3D variance for high sonic Mach number {M > 10). 
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